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ABSTRACT 



To improve photometric precision by removing blending effect, a newly 
developed technique of difference image analysis (DIA) is adopted by several 
gravitational microlensing experiment groups. However, the principal problem 
of the DIA method is that, by its nature, it has difficulties in measuring the 
baseline ffux Fq of a source star, causing degeneracy problem in determining 
the lensing parameters of an event. Therefore, it is often believed that the DIA 
method is not as powerful as the classical method based on the PSF photometry 
in determining the Einstein time scales ts of events. 

In this paper, we demonstrate that the degeneracy problem in microlensing 
events detectable from the searches by using the DIA method will not be as 
serious as it is often worried about. This is because a substantial fraction 
of events will be high amplification events for which the deviations of the 
amplification curves constructed with the wrong baseline fluxes from their 
corresponding best-fit standard amplification curves will be considerable even 
for a small amount of the fractional baseline fiux deviation AFq/Fq. With 
a model luminosity function of source stars and under realistic observational 
conditions, we find that ~ 30% of detectable Galactic bulge events are expected 
to have high amplifications and their baseline fiuxes can be determined with 
uncertainties AFq/Fq < 0.5. 
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1. Introduction 

Detection of a large number of events is one of the big challenges in microlensing 
searches. Classical solution to this challenge is observing fields with greatest density of stars 
such as the Galactic bulge (Alcock et al. 1997a; Udalski et al. 1997; Alard & Guibert 1997) 
and the Magellanic Clouds (Alcock et al. 1997b, 1997c; Ansari et al. 1996). While the use 
of such crowded fields increases the event rate, it also limits the precision of the photometry 
due to blending (Di Stefano & Esin 1995; Wozniak & Paczyhski 1997; Han 1997; Alard 
1997). In addition, with the use of the classical method based on PSF photometry one can 
monitor only stars with resolved images, and thus the number of source stars is limited by 
crowding. 

These problems of the classical method of microlensing experiments can be resolved 
with the newly developed technique of difference image analysis (DIA, Alard 1998, 1999; 
Alard & Lupton 1998; Alcock et al. 1999a, 1999b; Melchior et al. 1998, 1999). Since the 
DIA method detects and measures the variation of source star flux by subtracting an 
observed image from a convolved and normalized reference image, one can measure light 
variations even for unresolved stars. By using the DIA method, one can not only improve 
the photometric precision by removing the effect of blending but also increase the number of 
detected events by including unresolved stars into monitoring sources. In addition, the DIA 
method allows one to overcome the restriction of conducting lensing experiments toward 
only resolved star fields and thus can extend our ability to probe extra-galactic MACHOs 
(Gould 1995, 1996; Han 1996; Han & Gould 1996; Crotts & Tomaney 1996; Tomaney & 
Crotts 1996; Ansari et al. 1997, 1999). 

However, the principal problem with the DIA method in microlensing experiments 
is that, by its very nature, it has difficulties in measuring the unamplified flux (baseline 
flux, Fq) of a source star. This is because the observed light curveQ of a microlensing event 
obtained by the DIA method, -FbiA, results from the combination of the true amplification 
Aq and the baseline flux, i.e. 



'DIA 



F-F,ef = Fo(Ao-l), (1.1) 



where F and Fj-^i represent the source star fluxes measured from the image obtained during 
the progress of the event and the reference image, respectively. One significant consequence 
of this problem is that it produces degeneracy in determining the lensing parameters of the 
event (see § 3) like the degeneracy problem for a blended event whose light curve results 
from the combination of Aq and the blended flux. Therefore, it is often believed that the 
DIA method is not as powerful as the classical method based on the PSF photometry in 
determining the Einstein time scale ^e of an event. 



^Throughout this paper, we use the term 'light curve' to designate the changes in the flux 
of a source star, while the term 'amplification curve' is used to represent the changes in the 
amplification of the source star flux. 
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In this paper, we demonstrate that the degeneracy problem in microlensing events 
detectable from the searches by using the DIA method will not be as serious as it is often 
worried about. This is because a substantial fraction of events will be high amplification 
events for which the deviations of the amplification curves constructed with the wrong 
baseline fluxes from their corresponding best-fit standard amplification curves will be 
considerable even for a small amount of the fractional baseline flux deviation AFq/Fq. With 
a model luminosity function of source stars and under realistic observational conditions, we 
find that ~ 30% of detectable Galactic bulge events are expected to have high amplifications 
and their baseline fluxes can be determined with uncertainties AFq/Fq < 0.5. 



2. Mis-normalized Amplification Curves 

The standard form of the amplification curve of a gravitational microlensing event is 
related to the lensing parameters by 



Ao{u) 



u^ + 2 



m(m2 + 4)V 



2' 



U 



Z b-f^ 



tE,0 



+ P'o 



1/2 



(2.1) 



where u is the lens-source separation normalized in units of the angular Einstein ring 
radius 9e, and the lensing parameters Po, tmax, and tE,o represent the impact parameter 
for the lens-source encounter, the time of maximum amplification, and the Einstein ring 
radius crossing time (Einstein time scale), respectively. Once these lensing parameters are 
determined from the amplification curve, one can obtain information about the lens because 
the Einstein time scale is related to the physical parameters of the lens by 
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Ve 
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D,. 



1/2 



(2.2) 



where te = DoiOe is the Einstein ring radius, v is lens-source transverse speed, M is the 
mass of the lens, and Dou Di^, and Dos are the separations between the observer, lens, and 
source star. 

However, if the baseline source star flux of an event is misestimated by an amount 
AFo, the resulting amplification curve A (hereafter 'mis- normalized' amplification curve) 
deviates from the true amplification curve Aq by 



A 



FqAo + AFo Ao + f 



Fo + AFo 



1 + / 



(2.3) 



where / = AFq/Fq is the fractional deviation in the determined baseline fiux.0 If the 
baseline flux is overestimated (i.e. / > 0), the determined amplification is lower than Aq, 



^We note that if / represents the blended light fraction, i.e. / = B/Fq, equation (2.3) 
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and vice versa. Note that while there is no upper hmit for /, it should be greater than —1 
(i.e. / > -1). 

The shape of a microlensing event amplification curve is characterized by its height 
(peak amplification) and the width (event duration), which are parameterized by the impact 
parameter and the Einstein time scale, respectively. Since both the height and width of 
the amplification curve are changed due to the wrong estimation of the baseline flux, the 
lensing parameters determined from the mis-normalized amplification curve will differ from 
the true values. First, the change in the peak amplification makes the determined impact 
parameter change into 



/3 



2 1 



,-.^-V^ 
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1/2 
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Ap,o + f 
1 + / ' 



(2.4) 



where ApQ = (/9q + 2)/Pq{Pq + 4)-*^/^ and Ap are the peak amplifications of the true and the 
mis-normalized amplification curves. In addition, due to the change in the event duration, 
the determined Einstein time scale differs from the value tE,o by 



^E — ^E,0 



PI 



PI 



1/2 



Pth,0 



P' 



(2.5) 



where f3th represents the maximum allowed impact parameter (threshold impact parameter) 
for a source star to be detected by having a peak amplification higher than a certain 
threshold minimum value At^. With the right choice of the baseline flux, the required 
minimum peak amplification and the corresponding maximum impact parameter are 
Athfi = 3/v^ and Pthfl = 1- However, since the detectability will be determined from the 
mis-normalized amplification curve, the actually applied threshold amplification and the 
corresponding impact parameter will differ from Athfl and (3th,o by 



A 



th 



A 



th,0 



;i + /)-/, 



and 



(Han 1999). 
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(2.6) 



(2.7) 



In the upper panels of Figure 1, we present four example mis- normalized amplification 
curves A (solid curves) which are expected when the baseline flux of the source star for a 



describes the observed amplification curve of a microlensing event affected by blended light 
of an amount B. Therefore, the amplification curve of a blended event can be regarded as the 
mis-normalized amplification curve constructed with the baseline flux deviation AFq = B. 
The only difference is that since the blended light should be positive, i.e. / > 0, while the 
baseline flux deviation can be either negative or positive, blended amplification curves are 
always underestimated. 
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microlensing event with (3q = 0.5 is determined with the fractional deviations of / = ±0.2 
and ±0.5. By using equations (2.4) - (2.7), we compute the lensing parameters of the 
standard amphfication curves which best fit the individual mis-normalized amplification 
curves, and the resulting amplification curves Agt are presented by dotted lines. In the 
lower panels, to better show the difference between each pair of curves A and Agt, we also 
present the fractional deviations of the amplification curves A from their corresponding 
best-fit standard amplification curves, i.e. AA/Af^i; A A = Agt — A. From the figure, one 
finds the following trends. First, for the same amount of |/| = |AFo| /Fq, the fractional 
deviation AA/Af^ is larger when the baseline flux is underestimated (i.e. / < 0) compared 
to the deviation when the baseline flux is overestimated (i.e. / > 0). Second, although the 
difference between the two amplification curves A and Agt becomes bigger as the deviation 
AFq/Fq increases, the mis-normalized amplification curves, in general, are well fit by 
standard amplification curves with different lensing parameters. 



3. Baseline Flux Determination for High Amplification Events 

In previous section, we showed that since the amplification curve of a general 
microlensing event obtained based on wrong estimation of the baseline flux is well fit by 
a standard amplification curve with different lensing parameters, making it difficult to 
determine Fq from the shape of the observed light curves. In this section, however, we show 
that for high amplification events the deviations of the mis-normalized amplification curves 
from their best-fit standard curves are considerable even for a small fractional deviation 
AFq/Fq, and thus one can determine the baseline fluxes with small uncertainties. 

To demonstrate this, in the upper panels of Figure 2, we present the mis-normalized 
amplification curves constructed with the same fractional baseline flux deviations of 
/ = ±0.2 and ±0.5 as the cases in Figure 1, and the corresponding best-fit standard 
amplification curves for a higher amplification event with an impact parameter of f3o = 0.1. 
In the lower panels, we also present the fractional differences AA/Af^. From the comparison 
of the fractional differences AA/Af^i in Figure 1 and 2, one finds that the deviations of the 
mis-normalized amplification curves from their corresponding standard amplification curves 
are siginificantly larger for the higher amplification event. 

To quantify how better one can determine the baseline flux with increasing event 
amplifications, we determine the uncertainty ranges of Fq for microlensing events with 
various impact parameters /3o under realistic observational conditions. To do this, for a 
given event with f3o we first produce a series of mis-normalized amplification curves which 
are constructed with varying values of /. In the next step, we obtain the best-fit standard 
amplification curves corresponding to the individual mis-normalized amplification curves by 
using the relations in equations (2.4) - (2.7). We then statistically compare each pair of the 
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amplification curves A and Agt by computing x^, which are determined by 



x^=E 



1=1 



pAfitiU 



(3.i: 



For the computation of x^j we assume that the events are observed A^dat = 60 times 
during —Its < t < Its- The photometric uncertainty p depends on the observational 
strategy, instrument, and source star brightness. Therefore, we determine the photometric 
uncertainty by computing the signal-to-noise ratio (1/p = S/N) under the assumption that 
events are observed with a mean exposure time of tgxp = 150 s by using a 1-m telescope 
equipped with a detector that can detect 12 photons/s for a / = 12 star. The detailed 
description about the signal-to-noise computation is described in § 4. Once the values of x^ 
as a function of / are computed, the uncertainty of Fq is determined at la (i.e. x^ = 1) 
level. We then repeat the same procedures for events with various values of f3o (and thus 
the peak amplifications). 

In the upper of Figure 3, we present the resulting values of x^ as a function of log(l + /) 
for example events with the source star brightness / = 18 and various impact parameters 
of f3o = 0.05, 0.1, 0.2, and 0.22. In the lower panel, we also present the uncertainty range 
of AFq/Fq (shaded region). From the figure, one finds the following trends. First, the 
uncertainty significantly decreases as the impact parameter decreases, implying that the 
baseline fluxes for high amplification events can be determined with small uncertainties. 
Second, if the impact parameter becomes bigger than a certain critical value (/^crit), the 
value of x^ becomes less than 1, implying that the Fq cannot be determined from the 
shape of the obtained light curve. For our example events with / = 18, this corresponds to 
/^crit = 0.22. Note that the uncertainty range AFq/Fq in the lower panel is determined only 
for impact parameters yielding x^ ^ 1- Third, the upper limit of the uncertainty range is 
always bigger than the lower limit. 

Knowing that Fq can be determined only for high amplification events, we define the 
critical impact parameter f3cnt as the maximum allowed impact parameter below which the 
baseline flux of the event can be determined with uncertainty less than 50% (i.e. x^ ^ 1 
and AFq/Fq < 0.5. Then /3crit(-^o) represents the average probability that the baseline flux 
of an event with a source star brightness Fq can be determinded with an uncertainty less 
than 50%. We compute the critical impact parameters for events expected to be detected 
towards thd Galactic bulge, and they are presented in the upper panel of Figure 4 as a 
function of the source star brightness in / band. From the flgure, one flnds that as the 
source star becomes fainter, the value of /^cnt decreases. This is because for a faint source 
event, the photometric uncertainty p is large. Therefore, to be distinguished from standard 
ampliflcation curves with a statistical confldence level higher than the required level (i.e. 
X^ ^ 1)) the event should be highly amplifled. 



4. Fraction of High Amplification events 

In previous section, we showed that the baseline fluxes of high amplification events 
can be determined with precision. In this section, we determine the fraction of high 
amplification events for which one can determine Fq with small uncertainties among the 
total microlensing events detectable by using the DIA method. 

Under the assumption that image subtraction is perfectly conducted^, the signal 
measured from the subtracted image by using the DIA technique is proportional to the 
variation of the source star fiux, i.e. 5* oc Fq{Aq — l)texp- On the other hand, the noise of the 
source star fiux measurements comes from both the lensed source star and blended stars, 
i.e. N oc (FqAq + BY^'^, where B is the background fiux from blended stars in the effective 
seeing disk (i.e. the undistinguishable separation between images) with a size (i.e. diameter) 
A6'see. Then the signal-to-noise ratio of an event whose light variation is detected by using 
the DIA method is given by 

where (B) represents the mean background flux. For a high ampliflcation event {Aq ^ 1) 
with a bright source star {FqAq ^ (B)), the signal-to-noise ratio becomes photon limited, 
i.e. S/N ~ (Fo^otexp)^''^. By contrast, for a low ampliflcation event with a faint source 
star {FqAq <C (-B)), the noise from the background flux becomes important. Let us deflne 
/9max(-fo) as the maximum impact parameter within which a lensing event can be detected 
by having signal-to-noise ratios higher than a certain threshold value (S'/A^)th during a range 
of time longer than a required one At. Then, /9max(-Po) represents the average detection 
probability of an event with a source star brightness Fq from the microlensing search by 
using the DIA method, and it is computed by 

/Smax = I [""^^^ ~ i^t/^E?]^^^ when W„,ax > At/te (42) 

1 when u^ax < At/ts 

where Wmax = [2(1 — A:^^^)^^^"^ — 2]^/^ represents the threshold lens-source separation below 
which the signal-to-noise ratio of an event becomes greater than {S/N)th and Ay^in is the 
ampliflcation of the event then M=Mmax- The value of Amin is obtained by numerically 
solving equation (4.1) with respect to the ampliflcation for a given threshold signal-to-noise 
ratio of {S/N)t^. 

In the upper panel of Figure 4, we present the maximum impact parameter /3max as 
a function of the source brightenss for stars in the Galactic bulge. For the computation 



^ A very ingeneous image subtraction method developed by Alard & Lupton (1998) 
demostrate that it is possible to measure the variable flux to a precision very close to the 
photon noise. 
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of PraaxiFo), we assume the same observational conditions described in § 3. The adopted 
threshold signal-to-noise ratio of the MACHO experiment is (S'/A^)th = 10 (Alcock 1999a, 
1999b). In our computation, however, a higher value of (S'/A^)th = 15 is adopted to account 
for the additional noise from the sky brightness and the residual flux due to imperfect image 
subtraction. The average background flux is obtained by 

{B) = / FoMFo)dFo, (4.3) 

Jo 

where Fcl and $o(-fo) are the crowding limit of the Galactic bulge fleld and the luminosity 
function of stars in the fleld normalized to the area of 7r(A6'see/2)^. We adopt the luminosity 
function of Holtzman et al. (1998) constructed from the observations of bulge stars by 
using the Hubble Space Telescope and the adopted crowding limit is / = 18.2 mag. We 
assume that an event is detectable if signal-to-noise ratios are higher than (S'/A^)th during 
At = 0.2tE of the source star flux variation measurements. From the flgure, one flnds that 
the detection probabilities (i.e. /Smax) of events with source stars faint than the crowding 
limit, and thus unresolvable, are not negligible up to ~ 3 mag below Fch implying that a 
substantial fraction of events detectable by using the DIA method will be faint source star 
events (Jeong, Park, & Han 1999). 

With the determined values of /^max as a function of source star brightness, we then 
construct the effective source star luminosity function by 

$eff(Fo)=/5max(Fo)*o(i^o). (4.4) 

We also construct the luminosity function of source stars for high ampliflcation events with 
measurable baseline fluxes by 

dS — r,( rr \(h ( rr \- r, — ) Pcmt/Pmux when Pcnt ^ Pmax (A c-\ 

$high-^(Fo)$o(Fo), V-[^^ when /3e.t > /5.ax ^^-^^ 

Once the two luminosity functions $cflf and $high are constructed, the fraction of high 
amplification events out of the total number of detactable events is computed by 

^high ^ /o°°'^high(-^o)c^-^0 / , „N 

Ttot /o"$eff(i^o)rfi^O ■ ^'^ 

We find that for ~ 33% of events detectable from the microlensing searches by using 
the DIA method will have high amplification events for which the baseline fiuxes can be 
determined with uncertainties AFq/Fq < 0.5. 



5. Conclusion 

We have investigated how the lensing parameters change due to the wrong determination 
of the baseline fiux of a microlensing event. We have also investigated the feasibility of the 
baseline fiux determination from the shape of the observed light curve. The results of these 
investigations are as follows: 
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1. The obtained amplification curve of a general microlensing event based on wrong 
baseline flux is well fit by a standard amplification curve with different lensing 
parameters, implying that precise determination of Fq from the shape of the observed 
light curve will be difficult. 

2. However, for a high amplification event, the mis-normalized amplification curve 
deviates from the standard form by a considerable amount even for a small fractional 
deviation of baseline flux, allowing one to determine Fq with a small uncertainty. 

3. With a model luminosity function of Galactic bulge stars and under realistic 
observational conditions of the microlensing searches with the DIA method, we find 
that a substantial fraction (~ 33%) of microlensing events detectable by using the 
DIA method will be high amplification events, for which the baseline fluxes of source 
stars can be determined with uncertainties AFq/Fq < 50%. 



This work was supported by the grant (1999-2-113-001-5) of the Korea Science & 
Engineering Foundation. 
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Figure 1: Upper panels: four example mis- normalized amplification curves (solid curves) 
which are expected when the baseline flux of the source star for a microlensing event 
with Pq = 0.5 is determined with the fractional deviations of / = ±0.2 and ±0.5. Also 
presented are the standard amplification curves (Agt, dotted curves) which best fit the 
individual mis- normalized amplification curves. Lower panels: the fractional deviations 
of the mis-normalized amplification curves from their corresponding best-fit standard 
amplification curves, i.e. AA/Af[t', AA = Agt — A. 
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Figure 2: Upper panel: the mis-normalized amplification curves with the same fractional 
baseline flux deviations of / = ±0.2 and ±0.5 as the cases in Figure 1, and the corresponding 
best-fit standard amplification curves for a /i«5^/ier amplification event with /5o = 0.1. Lower 
panel: the fractional differences between the amplification curves A and Afif 
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Figure 3: Upper panel: the values of x^ as a function of log(l + AFq/Fq) for example 
events with a source star brightness of / = 18 and various impact parameters /5o- The value 
of x^ is computed by comparing the mis-normalized amplification curve with a fractional 
baseline flux deviation / = AFq/Fq and the corresponding best-flt standard ampliflcation 
curve. Lower panel: the uncertainty range of the baseline flux (shaded region) determined 
from the shape of the light curve of a lensing event detected by using the DIA technique. 
The uncertainties are determined at Icr (i.e. x^ = 1) level. 



-15- 





i 


CD 




*-> 

a; 


0,8 


e 




to 


0.6 


to 




Oh 


0.4 


^J 




o 




to 


?, 


Uh 




H 









CO 


800 



R 



o 

CO 



w 

CO 



C 
o 

•r— I 

o 

c 

-I — I 

m 
O 

CI 



600 



400 



200 




Figure 4: Upper panel: the critical and the maximum impact parameters (/3crit and /5max) 
as functions of the source brightness for stars in the Galactic bulge field. The value of /5max 
is equivalent to the average detection probability of an event with a source star brightness I 
from the microlensing search by using the DIA method. On the other hand, j3crit represents 
the average probability that the baseline flux of an event with a source star brightness / 
can be determined with an uncertainty less than 50%. Lower panel: the effective source 
star luminosity functions of the total ($eff) and high amplification events ($high) detectable 
from the microlensing searches by using the DIA technique toward the Galactic bulge field. 
The shaded region represents the fraction of high amplification events. 



